Coherent Collisions between Bose-Einstein Condensates 
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We study the non-degenerate parametric amplifier for matter waves, implemented by colliding 
two Bose-Einstein condensates. The coherence of the amplified waves is shown by observing high 
contrast interference with a reference wave and by reversing the amplification process. Since our 
experiments also place limits on all known sources of decoherence, we infer that relative number 
squeezing is most likely present between the amplified modes. Finally, we suggest that reversal 
of the amplification process may be used to detect relative number squeezing without requiring 
single-particle detection. 
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Quantum entanglement has always attracted great in- 
terest as it is one of the key differences between quantum 
and classical physics. Since it can be exploited for various 
novel applications like quantum computation and preci- 
sion measurements, there has been a continuing effort 
to engineer ever more robust quantum correlated states. 
Recently, four-wave mixing using Bosc-Einstcin conden- 
sates (BEC's) has been proposed as a method for creating 
two correlated matter waves fit]. In that experiment, two 
source waves and a third seed wave were used to pro- 
duce a fourth wave in the conjugate momentum state, 
resulting in an amplified seed-conjugate wave pair. For 
large amplification, the relative number fluctuations be- 
tween the amplified waves can become significantly re- 
duced, and this was recently achieved S. In optics, this 
is known as non-degenerate parametric amplification |3j , 
where it is also used for generating similarly squeezed 
photon beams. Once created, such two-mode number 
squeezed states could be used for Heisenberg-limited in- 
terferometry where phase shifts up to 1 /N accuracy could 
be measured [0, |j| . As the occupation number N of the 
states increases, the advantage gained over the standard 
shot noise limit of 1/VN becomes quite significant. 

Entanglement between the generated matter waves is 
possible only if the collisional four-wave mixing process 
is coherent. While the coherence of collisions has been 
observed before in mean field effects [[| ^ , it has not been 
demonstrated for the more common case of momentum 
changing elastic collisions. By itself, four wave mixing 
does not require or guarantee coherence, since it will still 
occur for condensates with short coherence lengths ||. 
Moreover, even though previous four wave mixing ex- 
periments |jL § utilized two coherent condensates, they 
did not explicitly show that coherence was maintained 
throughout the collisional process. 

In this paper, we demonstrate the coherence of elas- 
tic collisions by observing high contrast interference be- 
tween the conjugate wave and a reference wave. Due to 
the strong coherence retained, we also observe a reversal 
of the amplification process. The success of these two 
experiments places an upper bound on the decoherence 



effects in our experiments, providing evidence that a sig- 
nificant amount of relative number squeezing is present 
in the pair-correlated beams. The reversibility of the 
amplification process further suggests a possible way of 
directly detecting the two-mode number squeezing and 
entanglement between the two correlated states, without 
requiring sub-shotnoise detection. 

In the Heisenberg picture, the evolution of the annihi- 
lation operators a 3l of a mode pair that is parametri- 
cally amplified by the source modes ax, ai is given by 
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where V is the volume, m is the mass, a is the s-wave 
scattering length and Ae is the energy mismatch between 
the source modes and the amplified modes M. The 
solution to these equations for highly occupied source 
modes is an exponential growth with a rate constant 

7/ = yj (^) 2 — (^f ) 2 where ~p is the geometric mean of 
the chemical potentials of the two source waves. In our 
experiment, the multiple waves used were generated from 
a single BEC using two-photon stimulated optical tran- 
sitions (9). The evolution of annihilation operators a, b 
for modes which are coupled by such a Bragg transition 
is given by 
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where is the coupling strength. In both the Bragg and 
amplification process, the generated wave incurs an over- 
all phase shift of 90° relative to its source, a detail which 
becomes important for the interference experiment. 

The experimental setup was almost identical to the 
one used in Sodium condensates of about 13 mil- 
lion atoms each were created in a cylindrical magnetic 
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trap with an axial frequency of 20 Hz and a radial fre- 
quency of 40 Hz. They had a Thomas- Fermi radius of 
30 firm radially and 60 fim. axially, a chemical potential 
of h x 1.8 kHz and a speed of sound of 6 mm/s. From 
an initial condensate ki , a small seed k3 with a velocity 
of 14 mm/s was generated using a 20 [is Bragg (seed) 
pulse (see Fig. |l|). In the next 40 fis, a second source 
wave k2 with a velocity of 20 mm/s was split off from 
k2 using a it/ 2 pulse. All Bragg beams in our experi- 
ment were detuned from the sodium D2 line by 40 GHz 
to avoid spontaneous Rayleigh scattering. Subsequent 
collisions between ki and k2 acted as a parametric am- 
plifier for the seed wave k3, creating a fourth conjugate 
wave k4 in the momentum conserving state. Both k3 and 
k4 modes then grew exponentially with an expected rate 
of 17 ms _1 on average accross the condensate. At the 
same time, the parametric amplifier also enhanced spon- 
taneous scattering into initially unseeded modes, giving 
rise to the collisional halo seen in Fig. |] |l(J. The result- 
ing loss of atoms from the source waves limited the four 
wave mixing process to 400 /tzs. 

After 80 /is of amplification, we probed the phase co- 
herence of the conjugate wave k4 by interfering it with 
a reference wave. For optimal interference, the reference 
wave needed to have a well defined relative phase and 
be identical in size and momentum to k4. The refer- 
ence wave was outcoupled from k 2 using another 20 //s 
Bragg (seed) pulse. By changing the phase of one of the 
Bragg beams, the phase of the reference wave relative 
to k4 was varied. The resulting occupation of the con- 
jugate momentum state k4 was monitored by selectively 
transferring a fixed fraction of atoms in k4 to a detection 
state k' 4 (see Fig. |l|). For this, a 40 ^is Bragg pulse with 
large momentum transfer (2 hk) was used. This ensured 
that the k 4 atoms were well-separated during ballistic ex- 
pansion and their number was easily determined. Direct 
number measurement of the mode k4 was not possible 
since the amplification process continued during ballistic 
expansion. 

The interference fringe obtained is plotted in Fig. 0. By 
fitting a sinusoidal curve to the raw data, we obtained a 
fringe contrast of 68 ± 11 %. This proves that the conju- 
gate wave maintains the phase relation established during 
amplification and is in a single quantum state. In addi- 
tion, we expected destructive interference to take place 
at 90° due to the accumulation of all the 90° phase shifts 
introduced by the four wave mixing and Bragg processes 
(see eqs. |l|, ||). Our value of 77 ± 9° agrees fairly well 
with this, and the small discrepancy could be due to the 
Bragg beams not being fully resonant. It is also inter- 
esting to note that the position of the fringe minimum 
is a direct signature of the positive sign of the scattering 
length a, since a negative scattering length would have 
resulted in a minimum at 270°. Usually, only the cross 
section, which is proportional to a 2 , is measured using 
these momentum changing collisions. 




FIG. 1: The two source waves ki, k2, the seed k3, the conju- 
gate wave k4, and the detection state k 4 after 20 ms ballistic 
expansion. Also seen is the collisional halo around all four 
waves as well as several small unlabelled waves, generated by 
off-resonant scattering. The field of view is 1.0 x 1.9 mm. 
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FIG. 2: Coherence of the conjugate wave. Shown is the in- 
terference of the conjugate wave with a reference wave. The 
relative phase shift was introduced by shifting the phase of 
the rf-wave which drives the lower-frequency acousto-optic 
modulator in the Bragg setup. 



A direct implication of the sustained coherence is the 
reversibility of the amplification. If Ae is small, the sec- 
ond term in eq. |l| becomes negligible and the overall evo- 
lution is solely determined by the unitary four- wave mix- 
ing term. By changing the sign of one of the four waves, 
we can apply the inverse operator and the amplified pop- 
ulations will subsequently deamplify to zero due to the 
entanglement between the coupled modes. This is true 
even for the modes which were not deliberately seeded, 
since spontaneous emission into these modes also occurs 
by the same pairwise collision process. A similar experi- 
ment using photons has already been performed []lT| . 

After 80 (is of amplification, we initiated deamplifi- 
cation by inducing a 180° phase shift in k2. The sign 
reversal was achieved by applying a 120 /is 7r-pulse reso- 
nant with the momentum difference between the source 
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FIG. 3: Phase coherent amplification and deamplification: (•) 
Deamplification of the conjugate wave, (o) Continued ampli- 
fication due to a 7r-pulse with a 180° phase shift applied in 
mid-pulse (see text), (x) Uninterrupted amplification. The 
solid, dashed and dotted lines are independent numerical sim- 
ulations of the respective processes. 



waves ki and k2. (For this experiment, also the 7r/2-pulse 
which originally generated k2 was prolonged to 60 /is.) 
To check that the amplification process was not perma- 
nently disrupted by the Bragg pulse, we also performed 
a control experiment where the 7r-pulse was replaced by 
a 7r/2 and — 7r/2 pulse applied in quick succession. Like 
the first half of the 7r-pulse, the 7r/2 pulse transferred 
all the atoms from k2 back into ki, but unlike the sec- 
ond half of the 7r-pulse, the —ir/2 pulse built up mode 
occupation in k2 with the original phase. If the Bragg 
pulses did not otherwise disturb the phase relation of the 
system, amplification should then resume. In all three 
cases, we monitored the evolution of by applying the 
detection pulse and counting the number of atoms in k^ 
after variable intervals. 

Fig. ^ shows the results obtained and clearly demon- 
strates the reversal of the amplification process. As a 
quantitive measure of the deamplification process, we de- 
fined the deamplification factor as the ratio of the maxi- 
mum mode occupation to the remnant left at the end of 
deamplification, and obtained a factor of ~ 5. In addi- 
tion, we infer that this reversal is phase-dependent, since 
amplification continues if we do not reverse the phase 
of the second source wave (Fig. ||(o, x)). This demon- 
strates that the four-wave mixing process does not suffer 
significant decoherence during the observation time. 

In Fig. ||, we also plot numerical simulations of all three 
cases and find good qualitative agreement with the ex- 
periment. For these simulations, the energy mismatch 
Ae was the main adjustable parameter and a\ and 02 
were taken to be c-numbers since they are highly occu- 
pied relative to the other modes. Apart from Ae, the only 
other adjustable parameters in the simulations were an 
overall normalization constant and the size of the initial 
seed, which was assumed to be small. Numerically, opti- 



mal deamplification was achieved for Ae = h x 1.0 kHz, 
and adjustments to this value by h x 0.6 kHz either way 
still yielded a fit consistent with the experimental data. 
The simulations also took into account the off-resonant 
scattering of the seed and conjugate wave during the ir- 
pulse, which caused a 20% loss of atoms from the ampli- 
fied waves and accounted for the dip in mode occupation 
seen in Fig. 3(o). Finally, beyond 300 fis, amplification 
slows down in Fig. |3|(°) and departs from the numerical 
simulation (dashed line) due to depletion of the source 
waves as a result of scattering and subsequent amplifica- 
tion of the unseeded modes || . 

Theoretically, maximum deamplification should occur 
for Ae = in the absence of any extraneous phase shifts. 
Our experiment deviates from this ideal situation as it 
suffered from off-resonant scattering during the 7r-pulse 
resulting in a calculated phase shift of 0.6 rad in the 
amplified modes. This was due to the finite duration of 
the 7r-pulse, which off-resonantly coupled out atoms from 
the amplified waves. The observed value of Ae = h x 
1.0 kHz therefore might have compensated for this phase 
shift as the amplified modes evolved over time. 

At the same time, we also discovered that the deam- 
plification process is extremely phase-sensitive. Com- 
pared to the tolerance in Ae for amplification [Q, the 
h x 0.6 kHz tolerance for deamplification is relatively 
small. While amplification is inherently a stable process 
which enforces a certain phase relationship between the 
four waves, deamplification is unstable with respect to 
the phase. As soon as the phase relation is lost, deam- 
plification ends and amplification resumes. Fortunately, 
deamplification remains possible since its tolerance is still 
larger than the estimated h x 0.4 kHz halfwidth of the 
inhomogeneous mean field shift. However, in order to 
remain within the narrow range of Ae allowed, we had 
to align the Bragg beams to within milliradians of the 
optimal value. 

In principle, the pairwise collisions of atoms during 
deamplification could be used to demonstrate relative 
number squeezing and entanglement between the two am- 
plified modes. If decoherence is absent, the mode occupa- 
tions after deamplification will reflect the relative number 
uncertainty in the initial waves. For coherent or classical 
waves, the remnant will follow Poissonian statistics and 
be of order y/~N, where N is the initial mode occupation. 
In contrast, for perfect two mode number squeezing, the 
remnant will always be exactly zero. Moreover, complete 
deamplification of the seed and conjugate waves can only 
occur if the sum of their individual phases remain well- 
defined throughout the process. These two conditions 
cannot be met for uncntanglcd waves, since their individ- 
ual number and phase are canonically conjugate variables 
and cannot simultaneously be well-defined. 

To illustrate this, we simulated the evolution of the 
conjugate wave for our experimental conditions, this time 
without an initial seed (Fig. ^a). Here, complete deam- 
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FIG. 4: Deamplification of perfectly squeezed modes with 
off-resonant scattering (solid line) . Identical process but with 
coherent waves (dashed line). The vertical scale is the expec- 
tation value of the mode occupation. 



plification does not occur as the squeezing has been de- 
graded by off-resonant scattering. However, it still dips 
below the shot noise level given by the minimum of the 
control graph (Fig. |]b) which assumes coherent waves. 
The control simulation is started at 160 /xs, immedi- 
ately preceding the 7r-pulse, using non-entangled coher- 
ent waves with the same unity mode occupation. To de- 
tect relative number squeezing, we need to obtain deam- 
plification factors of greater than y/N. 

Experimentally, direct detection of relative number 
squeezing in our strongly occupied beams is not possi- 
ble since it would require a much larger deamplification 
factor than what we currently observe. In this, we are 
limited by the inhomogeneous mean field and the deple- 
tion of the source waves. However, for initially unseeded 
modes, a deamplification factor of five is already evidence 
for sub-shot noise correlations (Fig. [|). This shows that 
by monitoring a collection of modes with a similar energy 
mismatch, detection of the entanglement is possible with 
current experimental parameters. Practically, however, 
the occupation in these modes was masked by surround- 
ing modes with a different energy mismatch. Further 
complicating detection, the momenta of the states were 
blurred out during the motion of the outcoupled atoms 
in the gradient of the mean field of the condensate. 

Finally, we briefly consider other possible sources of de- 
coherence in our amplifier and discuss how their effects 
may be minimized. Collisions between the source and 
amplified waves can compromise squeezing as atoms are 
scattered into or out of the amplified waves in a random 
way. As shown in Q, at low density, no stimulated scat- 
tering occurs and the spontaneous scattering rate was 
estimated to be significant only after the source waves 
are about to be depleted. By further manipulating the 
size of the condensate such that the amplified waves sep- 
arate out before the source waves are depleted, we can 



ensure that the squeezing will be preserved. 

A further study of the non-degenerate matter wave am- 
plifier should explore the issue of the "beam quality" of 
the squeezed waves. The phase errors incurred due to the 
inhomogeneous mean field induces spatially dependent 
phase shifts in the wavefunction of the squeezed waves, 
rendering them hard to isolate. Other sources of phase 
errors are not a concern as they appear to have little ef- 
fect on our experiment. In particular, phase fluctuations 
across the condensate (first observed by ^) are absent 
for our parameters [[l2[ |l3). Moreover, since the atoms 
travel at a maximum velocity of 20 mm/s for a duration 
of 300 /is, the required coherence length of 6 fiia is easily 
satisfied. 

The collisional matter wave amplifier studied here com- 
plements the amplifier based on optical superradiancc 
]l4] , While the gain and interference contrast of 

both are comparable, the collisional amplifier has a much 
larger momentum bandwidth, amounting to a significant 
fraction of the speed of sound In contrast, the optical 
amplifier has a bandwidth of only a single mode fl(| . 

In conclusion, we have demonstrated the coherence of a 
non-degenerate matter wave amplifier by measuring the 
phase of the conjugate wave and by reversing the am- 
plification process. While we do not yet directly detect 
number-squeezing in our amplified beams, we have evalu- 
ated all known decoherence processes and concluded that 
they are well controlled in our experiment. Therefore the 
atomic non-degenerate parametric amplifier studied here 
is a promising candidate for generating entangled and 
number squeezed beams. 
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